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In ultrathin ferromagnets deposited on metallic substrates, excitation of precessional motion of 
the spins produces a spin current in the substrate that transports angular momentum out of the film. 
This phenomenon is referred to as spin pumping, and is a source of damping of the spin motion. Spin 
pumping enters importantly in the description of spin dynamics in other nanoscale and subnanoscale 
systems as well. In this paper, we present an approach based on the Kubo formalism that allows the 
explicit calculation of this spin current and its spatial variation. We use the formalism to explore 
features of the spin current generated by spin motions in a simple model system. 

PACS numbers: 



I. INTRODUCTION 

Spin dynamics in magnetic nanostructures has been 
under active study for the past two decades. A very 
important issue is the nature of mechanisms that damp 
the spin motion. A central question is whether one en- 
counters new processes unique to the nanoscale environ- 
ment not encountered in bulk matter, where now damp- 
ing mechanisms are well understood. 

The answer to the question just posed is in the affirma- 
tive. An extrinsic mechanism referred to as two magnon 
scattering is evident in magnetic nanostrutures[T]. A sec- 
ond intrinsic mechanism, which is the topic of the present 
paper, is referred to as spin pumping. For this process to 
be active the spins of interest must reside in, or be cou- 
pled to a metallic structure. The most commonly studied 
system wherein spin pumping is observed is an ultrathin 
ferromagnetic film fabricated from a 3d magnetic metal 
or alloy that is deposited on a metallic substrate. Preces- 
sion of the spins within the ferromagnetic film generates 
a spin current in the substrate that transports transverse 
angular momentum out of the film, in the direction nor- 
mal to the interface between the film and the substrate. 
A consequence is that the amplitude of the precessional 
motion within the film decays with time. 

A theoretical discussion of the consequences of an- 
gular momentum transfer between precessing spins and 
conduction electrons was presented by Slonczewski[2], 
and Berger[3j explored relaxation of precessing moments 
by the spin pumping process. In an elegant discus- 
sion, Simanek and Heinrich[3] derived a correction to the 
damping term of the Landau-Lifshitz-Gilbert equation 
with origin in spin pumping. A short time later one of 
the present authors[5] and Simanek[Hj provided expanded 
discussions of spin pumping through use of the approach 
used in ref. 01 It should be remarked that the theoretical 
treatments just cited all employ a simple model wherein 
the precessing moments are described as localized spins, 
while the band electrons are treated as free electrons in 
a parabolic band. These theories thus outline the funda- 
mental physics associated with spin pumping, but their 
application to real materials is limited. The phenomenon 



of spin pumping was discovered experimentally by Urban, 
Woltersdorf and Heinrich[7]. 

In substrates wherein the spin diffusion length is very 
long, the spin current penetrates deeply. Thus, it can 
enter a second ferromagnetic film placed downstream to 
excite its spins through the spin torque mechanism[H]. 
One may envision devices made which exploit this phe- 
nomenon, most particularly when spins in the primary 
film are excited by short pulses. In a recent paper, the- 
oretical studies of pulse excitation of a second magnetic 
moment through spin current generated by a primary 
moment have been presented [9], in a closely related but 
somewhat different physical context than ultrathin ferro- 
magnets. 

While the theoretical discussions cited above clearly 
outline the basic physical origin of the spin pumping phe- 
nomenon, as noted they employ a very simple model of 
the band electrons so their application to real materials 
is limited. Within the framework of a semi-classical de- 
scription of spin motions in an ultrathin film sandwiched 
between two reservoirs, a description of spin pumping 
applicable to real materials has been developed [10] and 
calculations that account nicely for data have been pre- 
sented as well[TT]. In this approach, one electron theory 
forms the basis for the discussion. More recently a full 
quantum theoretic formulation of spin excitations in ul- 
trathin films in actual experimental geometries has been 
developed and applied to the analysis of ferromagnetic 
relaxation under circumstances where the spin pump- 
ing mechanism is the dominant contribution to the line 
widthjl2]. Here all features of the problem are described 
within itinerant electron theory, including the moments 
whose precession generates the spin current. The discus- 
sion is framed within the language of the Kubo formalism 
of many body theory. Successful accounts of spin pump- 
ing data were obtained in parameter free calculations for 
Fe on Au(100)[7], and also the theory accounts nicely 
for systematic features in the line width found in FMR 
studies of a FM/Cu/FM trilayer deposited on metallic 
substrates [15]. 

It is ferromagnetic resonance studies of spin waves in 
ultrathin films that have motivated our earlier theoret- 
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ical studies of the spin pumping mechanisms just cited. 
We point out that in ultrathin ferromagnets, spin polar- 
ized electron energy loss spectroscopy (SPEELS) can be 
used to excite and study large wave vector spin waves [H]. 
The line width of such modes is so very large that they 
exist only for two or three cycles of oscillation. The same 
spin pumping mechanism discussed in the literature on 
ferromagnetic resonance accounts quantitatively for the 
line widths found in SPEELS experiments, as discussed 
in ref. [12] 

It should be noted that none of the theoretical studies 
described above directly explore the spatial dependence 
of the spin current generated by the precession of the 
magnetization; as we shall see below, the underlying lat- 
tice has strong effects on its spatial structure. The ap- 
proach in ref. ITOland ref. [11] led to formulae only for the 
flux of spin angular momentum that flows out of the fer- 
romagnetic film, and in ref. [TO] it is the response function 
of the spin system that is calculated, with line width ex- 
tracted through a procedure very similar to that used by 
experimentalists when data is analyzed. The spin current 
itself is not calculated in this approach. 

It is important to develop the means to calculate the 
spin current directly, at the microscopic level and to ex- 
plore its spatial variation, including effects of the under- 
lying lattice. Model calculations presented below show 
these can be striking. It is noteworthy that tempo- 
ral dynamics and spatial distribution of pure spin cur- 
rents have been recently observed by nonlinear opti- 
cal effects[TOl [TO]. In this paper, we present a discus- 
sion within the framework of the Kubo formalism within 
which such questions can be explored. The approach may 
be regarded as extension of the methodology of ref. [12] 
and that used in our earlier studies of large wave vec- 
tor spin waves [T71 ITS] . All electrons in the system, in- 
cluding the moments whose precession is responsible for 
generation of the spin current, are described within the 
framework of itinerant electron theory appropriate to the 
3d based magnetic materials. Thus we do not resort to 
a picture within which the precessing spins that gener- 
ate the spin current are described phenomenologically. 
So in principle, we are not confined to the one electron 
description of the materials involved, in contrast to the 
discussions of ref. [10] and ref. [TT] We also present model 
calculations of the spatial dependence of the spin current 
in the presence of an underlying lattice, including the in- 
fluence of a surface on its behavior. This is done for a 
simple one-dimensional model within which we can ex- 
amine general aspects of the influence of the underlying 
lattice on the spin current. We remark that our cal- 
culations assume that the electrons propagate through 
the lattice in a ballistic manner, so our results apply 
to circumstances where the electron mean free path is 
long. Room temperature measurements[TO] of the indi- 
rect exchange coupling between Fe films separated by Au 
layers, all deposited by molecular beam epitaxy (MBE), 
show excellent agreement with theory based on ballistic 
electron propagation out to thirty lattice constants, be- 



yond which the exchange cannot be measured. In such 
high quality MBE grown samples, at low temperature 
one may expect ballistic descriptions of electron propa- 
gation, both in regard to the indirect exchange coupling 
discussed in ref. [TO] and the spin pumping currents ex- 
plored here, to extend to very long distances. We note 
also quantitative calculations of spin pumping contribu- 
tions to the linewidth that address real materials, such as 
those presented in refs. ITOHT2T employ ballistic descrip- 
tions of electron propagation through the lattice. 

In our view, interest in direct knowledge of the charac- 
ter of the spin current and its spatial distribution is stim- 
ulated by observations of the inverse spin Hall effect [2"U]. 
If coherent spin motions are excited in an ultrathin film 
on a metallic substrate, as discussed above, a spin current 
propagates normal to the interface between the film and 
the substrate. The spin current contains a DC compo- 
nent. If spin orbit coupling is present, then a DC charge 
current parallel to the interface is realized. This produces 
a Hall voltage across the film with magnitude propor- 
tional to the square of the amplitude of the spin preces- 
sion. The ability to perform explicit calculations of the 
spin current, and its spatial distribution is the first step 
in a microscopic theory of the inverse spin Hall effect. We 
note that descriptions of spin dynamics of nanomagnetic 
systems within which spin orbit coupling is present have 
appeared very recently [2~Tl I22j. In the future it should 
be possible to combine such discussions with analyses of 
spin currents such as described here to address the in- 
verse spin Hall effect. 

The outline of this paper is as follows. Section [ITJ is 
devoted to the theoretical formulation of how one gener- 
ates microscopic descriptions of spin currents within the 
framework of the Kubo formalism. Section Hill is devoted 
to presentation of a series of calculations for a model sys- 



tem, and concluding remarks are found in section IV 



II. THEORETICAL DESCRIPTION OF SPIN 
CURRENTS IN METALS 



A. General Comments 



It will be useful to provide a bit of orientation regard- 
ing the nature of spin currents in metals by recalling the 
principal results of the analysis found in ref. [TO] and ref. 
1111 These authors assume the magnetic moments in an 
ultrathin fcrromagnet precess coherently and in phase, so 
the magnetization is characterized by a rigidly rotating 
vector. We write this as M(t) = M s h(t) with M s the 
saturation magnetization and n(t) a unit vector [23 . The 
authors of ref. [TO] then show that the spin current radi- 
ated into the two semi infinite substrates which surround 
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a ferromagnetic film has the form 
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where we have taken h = 1. 

The expression in Eq. [I] involves vector quantities, 
whereas a complete description of the spin current re- 
quires a tensor object, with one index referring to the di- 
rection of transport and the second the spin polarization 
in the current. In the geometry considered, the direction 
of transport is normal to the surfaces of the film, and this 
is not noted explicitly. The vector character thus refers 
to the spin polarization carried by the current. The ex- 
pressions for A r and Ai will be discussed below. 

Suppose we consider excitation of the spins by a driv- 
ing field of frequency w, as in a ferromagnetic resonance 
experiment. Then in response to the field, we shall have 
M(t) = MsZ+5M(t), where in the linear response regime 
8M(t) oscillates with frequency to. We then have two 
components to the spin current. One is an AC compo- 
nent with the frequency uj 
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and then there is a DC component proportional in 
strength to the applied microwave power for small am- 
plitude spin motions. This is 



I 



j(DC) 



An M s 



dt 



(3) 



where the angular brackets denote a time average. It is 
the DC component that leads to the inverse spin Hall 
voltage discussed in the literature, when spin orbit cou- 
pling is present in the system. Notice, incidentally, that 
an AC Hall voltage also is necessarily present. 

In our analysis, within the framework of the Kubo for- 
malism, we generate an expression for the spin current. 
Our approach yields the AC component that, in the for- 
malism of ref. \W\ is described by the expression in Eq. 
[5J However, when an external driving field is applied 
to the moment or moments responsible for generating 
the spin current, we calculate Ig^^ directly without the 
need to resort to explicit descriptions of the one electron 
eigenstates of the system, and we calculate 5M(t) as well 
within the same theoretical structure. Thus, if desired 
we may extract the coefficients A r and Ai by comparison 
of our results with Eq. [2] which shows the polarization of 
the AC spin current is perpendicular to the static mag- 
netization. The coefficient Ai controls the magnitude of 
the component parallel to the time derivative of 5M{t) 
while A r controls the component perpendicular to this 



vector. With A r in hand, we may use Eq. [3] to generate 
an expression for the DC spin current. 

The point of these remarks is as follows. Our approach 
is based on linear response theory, and it generates the 
AC component of the spin current that is first order in the 
amplitude of the transverse magnetization, whereas the 
DC component is proportional to its square. Nonetheless 
if one accepts Eq. [T]as a description of the spin current, 
we can describe the DC component as well in our formal- 
ism. 



B. A Kubo Formalism Description of the Spin 
Current Generated by Application of an Exciting 
Microwave Field to a Moment Bearing Structure 

We now turn to a description of our method. In this pa- 
per, we confine our attention to a simple one dimensional 
lattice of sites, with one orbital per site. It is straightfor- 
ward to extend what follows to a three dimensional lat- 
tice with several orbitals per site by adding appropriate 
indices, but here we confine our attention to the simple 
system just described, in the interest of simplicity. 

The Hamiltonian that forms the basis of our treatment 
describes a lattice of sites occupied by atoms that need 
not be identical: 

H = ^ tijc\ a Cja- + ^ £ih la + i ^ Uifliahia 



(4) 
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One may view the Hamiltonian in Eq. [4] as describing a 
Hubbard alloy. In Eq. |4j tij is the hopping integral be- 
tween sites i and j, c\ a (ci a ) create (destroy) an electron 
with spin a at site i, n ia — ff ia Ci a , e.; is the energy of the 
orbital at site i and f/j represents the on site Coulomb 
interaction at site i. While our formal discussion allows 
Ui to be non-zero at every site in our system, in our nu- 
merical studies we shall set Ui to zero everywhere save 
for the site i — N where the magnetic moment is located 
whose precession generates the spin current we study. It 
is Ui which drives the formation of the magnetic moment 
at site N. Finally B is an externally applied magnetic 
field in the z direction, Sf = \{hi\ ~ fkl)- We refer to 
the magnetic moment at site N as too. 

We introduce the generalized spin operator = 

\ E b ^ia a a8^oP wnose diagonal component is the 
H th Cartesian component of spin in unit cell i. One finds 
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where If is the spin current radiated by precessing mag- 
netization in unit cell i. The contribution from the terms 
with k > i describe spin current propagating to the right, 
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and the contribution from terms with k < i describe spin 
current flowing to the left. If we consider a particular 
volume V that contains many spins, then the total spin 
current flowing out of V is 
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In Eq. [6] the sum over k may be restricted to k V 
because the contribution to l£ from the terms associated 
with k € V vanishes. It is useful to define the spin cur- 
rent operator Ig = Ig + ilg, from which one may easily 

obtain Ig and It. The amplitude of the total spin cur- 
rent pumped into the non-magnetic substrate due to the 
precession of rho is given by 



m)\ = \/m)\ 2 + m)\ 2 
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where (O) represents the expectation value of the oper- 
ator O. It is this quantity we study numerically in the 
discussion that follows. 

The applied transverse magnetic field b± (i) couples to 
the spin density of the magnetic unit located at site 0, 
and we shall calculate the change in the expectation value 
of the spin current caused by this relatively weak time de- 
pendent perturbation. Here, the change §(Ig) = 
because there is no net spin current flowing in the unper- 
turbed system. The interaction with a single frequency 
oscillatory b±(t) may be written as 



Hi = g^Bb±(t) ■ S = g/J, B bo cos(u)t)S% - sin(ut)S% 
g^Bbo ( &iut e _ iwt £_ 



(8) 

where the transverse field amplitude 6o is much smaller 
than B . 

Following the usual steps employed in linear response 
theory one finds that the spin current generated by this 
perturbation is given by 

<#>(*) = ^ E ^ l faxUM - , 



where 



(9) 
(10) 
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is the Fourier transform of the generalized susceptibility 

xtfkM = (0+W;^(o)» = -<e(t)([^(t),^(o)]) . 

(ii) 

Eq. [9] is the central result of this section. It relates the 
spin current pumped into the non-magnetic substrate to 
generalized transverse dynamical spin susceptibilities of 
the system, and we shall proceed next by briefly outlining 
how these susceptibilities may be determined. 



C. Generalized Transverse Spin Susceptibility 

By generating the equation of motion for XyjwW 
within the random phase approximation (RPA) one finds 

Xtjkli^) = Xijkli") ~ E Xij mm (u)U m Xmmkl(u) , (12) 
m 

where X%u ( a - 1 ) * s ^ ne noninteracting susceptibility, which 
may be written as 



v.*7ft«v-v - 2n I duJ f( u Hfr/U" ) ~ 9u ( w )]5,fc( w + w ) 
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(13) 



Here, gfj(ui) and g~° (lu) represent the time Fourier trans- 
forms of the retarded and advanced single-particle prop- 
agators for an electron with spin a between sites i and 
j, respectively, and f(ui) is the Fermi-Dirac distribution 
function. Direct numerical evaluation of the u/ integrals 



in Eq. 13 along the real axis may be difficult to be re- 
alized accurately because the one-electron propagators 
g and g~ usually are not smooth functions of to'. It is 
possible, however, to extend some those integrals to the 
complex plane and avoid a great deal of numerical dif- 
ficulties by rewriting Eq. [13] as a sum of three terms 
(x = h + h + h) which are given at zero temperature 

by m- 
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Once all the elements of x°( a; ) have been determined, 
Eq. [12] may be solved in matrix form and one may finally 
obtain the interacting susceptibility 

x + -H = [i + x°Ht>]- 1 x°H , (17) 

where U mn = U m 8 mn . 



III. NUMERICAL STUDIES OF THE SPIN 
PUMPING CURRENT 

In this section, we present the results of our studies 
of the spatial variation of the spin pumping current. We 
shall also compare the results obtained with our approach 
with those generated by the method described in ref. [10] 
and ref. 1111 We have used the simple one dimensional 
model from section In] for this first set of calculations. 
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FIG. 1: Modulus of the expectation value of the spin cur- 
rent }| produced by the precession of the moment at site 
0, calculated as a function of the energy E — fku across the 
boundaries shown in the inset. The spin current is in units 
of g(J,Bbo and the Zeeman field is adjusted so the resonance 
frequency of the impurity is Eo = hujo = 10 -3 . These calcu- 
lations are for the case where the band is half filled. 



We take nearest neighbor hopping only into account, 
and we set the hopping integral t = 1, so the bandwidth 
is 4 units. The zero of the energy scale is the middle of 
the band, so if the band is half filled the Fermi energy 
Ep = 0. We set the Coulomb interaction Ui to zero 
everywhere save for the magnetic moment bearing site 
where we take Uq — 10. We also set the orbital energies 
€i to zero everywhere save for site and eo is adjusted so 
the occupancy hq = 1. 

We begin by considering the moment embedded in a 
line of infinite length. The parameters mentioned above 
give a self-consistent moment on site of 0.916/jb. In 
this case, spin current is emitted by the precessing impu- 
rity both to the left and to the right. In Fig. [l] for the 
case where the band is half filled, we show the modulus 
of the expectation value of the total spin current, ) |, 
as a function of energy E — hu>, where uj is the frequency 
of the driving field, calculated across different bound- 
aries shown in the inset. The Zeeman field has been 
adjusted so that the resonance frequency of the impurity 
is Eo = hu>o = 10~ 3 . We see a resonance centered very 
close to Eo. Clearly there is spatial variation in the spin 
current, and also it is evident that the moment emits an 
appreciable spin current even for precession frequencies 
well above ujq. We remark that, as we see from Fig. [l] 
our Kubo formalism based approach allows us to directly 
calculate the absolute amplitude of the spin current as a 
function of the frequency that drives the spin system, 
with magnitude of the external driving field fixed. We 
do not need to link our analysis into the Landau-Lifshitz- 
Gilbert equation as in earlier studies [3HE] . 

In Fig. [2j we show the spatial variation of the spin 
current, for the case where the impurity is excited on 
resonance. In Fig. |2ja), the calculation is for the case 
where the band is half filled (Ep — 0). The magnitude 
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FIG. 2: We show the spatial dependence of the spin current 
for three choices of the Fermi energy, (a) Ef = (band half 
filled), (b) E F = 0.1, and (c) E F = 0.2. The impurity is 
driven on resonance in these calculations. 



of the spin current decays as one moves from the impu- 
rity, to approach a constant value at large distances. We 
see oscillations at half the Fermi wavelength about the 
average value of the spin current. When the band is half 
filled, the period of oscillations A — is twice the lattice 
constant. In Fig. [5Jb) and Fig. [2jc) we show the spatial 
variation for the case where the Fermi energy is moved 
off half filled. The Fermi wave vector for these two cases 
is incommensurate with the lattice, and we see the long 
wavelength oscillations that are very similar in nature 
to the phenomenon of aliasing that has entered earlier 
discussions of coupling between ferromagnetic films me- 
diated by a metallic non magnetic spacing layer [25]. We 
expect the phenomenon of "dynamic aliasing" illustrated 
m Fig. [2] will survive in calculations of the spin current 
generated by the coherent precession of spins in a film ad- 
sorbed on a substrate, since it is evident in calculations 
of the static interfilm couplings within the framework of 
three dimensional analyses [25]. 

It is of interest to compare the results of our calcula- 
tions with the picture set forth in ref. [inland ref. |TT| For 
the AC component of the spin current, Eq. [2] provides 
the relation 
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The amplitude of the spin motion may be generated 
from the Landau-Lifshitz-Gilbert equation that is modi- 
fied in form in the presence of spin pumping. The authors 
of ref. [TU] and ref. |TT] obtain the form 
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Here 



, , 7 and a 1 = 4 ""°"+g^ . 

7 = — gfig is the gyromagnetic ratio, mo is the mag- 
netic moment of the impurity in Bohr magnetons, and 
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a describes the intrinsic damping rate when spin pump- 
ing is absent. In our comparison between the two sets 
of calculations, we set a to zero. The magnetization 
is exposed to the magnetic field B(t) = B n z + b±(t), 
where b±(t) = bo [cos(wt)x — sin(wt)y]. One finds, after 
linearizing Eq. |19| 



-g^sbo 



Arf ~ \~ lAj 



iirmo + gAi 



TUqUJ 



la uj 



( 2 °) 

To proceed with the comparison, one needs to link the 
coefficients A r and A4 to our formalism. From ref. ITUl 
which employs a ballistic description of electron propa- 
gation as we do here, one has 

1 



. 1 . 



+ \t 



t 

mn 



ttnuV) (21a) 



and 



Ai Im ^ ^ { r mn{ r rnn) "f~ 4n(4n) } 



(21b) 



In these expressions, r^ nn and t^ nn are reflection and 
transmission coefficients for an electron with spin a at 
the Fermi energy Ep scattered by the magnetic unit 
from channel m to channel n of the normal metal. For 
our one dimensional model, there is only one channel 
available to both the initial state and the final state, so 
r mn — '■"tti^no an d similarly for t a mn . One finds [55] 
V ° 9 ™ {E ^ and t a = , T ,„ \ , where V a = 



l-V"g a 00 (E F ) 



l-V°g» (E F , 

-er^p with a = ±1 for majority (f) and minority 
(4.) spins, respectively, and <?o is the unperturbed one- 
electron propagator at site 0. 

In Fig. [3j we present a comparison between the results 
of our formalism, and the predictions that follow when 
Eq. [20] is used in conjunction with Eqs. 21 Both sets 
of calculations have been performed for identical model 
parameters. In Fig. [3]ja) we show the results generated 
with the formalism developed in the present paper. Here 
we evaluate the spin current very far away from the im- 
purity, where the spin current has settled down to its 
constant, asymptotic value. Fig. |3jb) gives the results 
from the formalism of ref. 10. It is the case that the 
calculations based on our formalism give a renormalized 
resonance frequency slightly higher than that which fol- 
lows from ref. [TUJ We suggest this has its origin in the 
fully dynamical nature of our analysis, wherein the spin 
current can react back onto the local moment and affect 
its precession frequency. We noted a similar effect earlier 
when full dynamical calculations of spin wave dispersion 
relations in ultrathin films are compared by those gen- 
erated by adiabatic theory [T7j. In the present case, the 
difference is rather small. 

If the precessing impurity is near a surface, the spa- 
tial variation of the spin current radiated from it differs 
dramatically from the case where it is in the bulk of an 
infinitely extended material. In the latter case, spin cur- 
rent streams both to the right and the left, and symme- 
try requires the amplitude of each stream to be the same. 
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FIG. 3: We show a comparison between (a) the spin current 
far from the precessing impurity calculated with our formal- 
ism, and (b) the spin current as described by the formalism of 
ref. fTOland ref. [11] The solid, dashed and dotted lines refer to 
calculations for which the Zeeman energy Eq — fiwo is taken 
to be 1.0 x 10~ 3 , 2.0 x 10" 3 and 3.0 x 10" 3 , respectively. 



However, as argued some time ago in a phenomenologi- 
cal discussion presented by Valet and Fert[27], the spin 
current must vanish at the surface. If the surface is lo- 
cated at z = L, then between our impurity and the sur- 
face, or in a system with a film radiating spin current 
into a substrate of finite thickness, in Valet-Fert theory 
the spatial profile of the spin current will have the form 



sinh ( 



where As is the spin diffusion length. In our 



model, the spin diffusion length is infinite, so within the 
phenomenological theory the spin current should have a 
linear variation, proportional to [L — z). 

An alternative explanation for the linear variation of 
the spin current is the following: at an arbitrary position 
z between the magnetic impurity and the surface, the 
total spin current I is given by the difference between 
the current I r flowing toward z — L and the current Ie 
flowing toward the impurity located at z = 0. Both I r 
and Ii are proportional to the direction of the transverse 
component of the electron spin at z, which is given by 
<x_l(t) = cos(wr)a; + sin(ojr)y, where r is the time the 
perturbation takes to reach z. Since the spin disturbance 
propagates with the electronic Fermi velocity, 77 = z/vf 
and r r = (2L — z)/vp. For relatively small frequencies, 
it is possible to show that |/| a \2u(z — L)/vp\. 

In Fig. [4j we show results of the spin current gener- 
ated by a precessing moment in a semi-infinite line of 
sites, calculated on resonance. The surface is located at 
site JV = 1, and the precessing impurity is located at 
site iV = 1001. We see that to the right of the precess- 
ing moment, the spin current settles down to a constant 
value as angular momentum is carried off into the bulk 
of the material. However, on the left side, in the region 
between the moment and the surface, our microscopic 
calculation produces a linear decrease of the spin current 
in excellent agreement with the phenomenology set forth 
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FIG. 4: The spatial variation of the spin current generated by 
a magnetic moment which precesses in a semi-infinite line of 
lattice sites, calculated on resonance. The surface is at N = 1, 
and the precessing impurity is located at the site N — 1001. 
In this calculation we set Ef — 0. 



in ref. 27. The surface influences the spatial variation 
in the spin current even when the precessing moment is 
very far from the surface. Of course, in a real material 
with finite spin diffusion length As the moment will not 
feel the presence of the surface when the distance be- 
tween it and the surface is long compared to As- In high 
purity materials such as the noble metals, the spin dif- 
fusion length can be quite long, it should be noted. One 
can grow ultrathin 3d ferromagnets on semiconducting 
substrates such as GaAs. It would be of great interest 
to cap such films with a noble metal such as Ag (indeed 
this is often done to prevent oxidation of the ferromag- 
net) and then study the dependence of the spin pumping 
relaxation rate as a function of thickness of the capping 
layer. 

From Fig. [4j in our model with its infinite spin dif- 
fusion length, it is clear that the presence of the surface 
strongly influences the nature of the spin current radi- 
ated from the precessing impurity. One can then inquire 
how, as the impurity is moved farther and farther from 
the surface, one approaches the spin pumping line width 
realized in the bulk of the material. This can be studied 
by generating the response function associated with the 
impurity moment itself. In principle, one can measure 
the width of the calculated resonant structure as a func- 
tion of distance of the moment from the surface, to see 
how the line width approaches the bulk value. 

We have studied this issue, to find that the approach 
to the bulk limit is not simple, for the following reason. 
When the impurity is at any given position N, the re- 
gion between site N and the surface at N = 1 acts as a 
quantum well of finite depth, by virtue of the fact that an 
electron reflected from the surface is backscattered by the 
one electron potential associated with the impurity site. 
There are then resonant levels for electrons in this region, 
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FIG. 5: The energy variation of the spin current radiated by 
an impurity located at position TV = 1001. The calculations 
show the influence of the quantum well states between the 
impurity position and the surface. Illustrated is the spin cur- 
rent calculated between sites 1001 and 1000 (solid line) and 
the spin current radiated off into the bulk, calculated between 
sites 1001 and 1002 (dashed line). The vertical dotted lines 
indicate the energies at which the results depicted in Fig. [6] 
are evaluated. 



and these produce dramatic structures in the resonance 
lineshape. Under these circumstances, it is difficult to de- 
fine the line width. We illustrate such structures in Fig. 
[5] where we plot the frequency variation of the spin cur- 
rent in the region between the impurity and the surface 
(solid line) and also the frequency variation of the spin 
current radiated into the bulk (dashed line), both cal- 
culated across the bond between the magnetic impurity 
and its nearest neighbors. Similar structures are found 
in the lineshape of the impurity response as well. As the 
impurity is moved farther and farther from the surface, 
these subsidiary peaks move inward to merge with the 
main structure in the resonant response. It is thus not 
clear how one can extract something that can be called 
a line width when the profile of the response function is 
so complex. It is the case that our calculation produces 
a bulk like response very far indeed from the surface, if 
one blurs out the fine structure produced by the resonant 
states. 

The spatial variation of the spin current in the region 
between the impurity and the surface is most striking, 
for excitation frequencies in the vicinity of the subsidiary 
resonances illustrated in Fig. [5] We illustrate this in Fig. 
[6] In Fig. [6]ja), we see that the spatial variation is con- 
trolled by a single resonance level (virtual level) in the 
region between the impurity and the surface, and the 
same is true in Fig. [6^b). A proper quantum theory is 
required for such features to appear. They are absent 
from the macroscopic approach described by Valet and 
Fert[27], for example. It will be of great interest to see 
if such features will be present in a full multiband de- 
scription of spin pumping relaxation in capped ultrathin 




FIG. 6: The spatial variation of the spin current generated 
by a magnetic moment which precesses in a semi-infinite line 
of lattice sites, calculated at energies (a) E — 0.0119 and (b) 
E = 0.0182, depicted by vertical dotted lines in Fig. [5] The 
surface is at N — 1, and the precessing impurity is located at 
the site N = 1001. In this calculation we set Ep = 0. 



FIG. 7: The spatial variation of the spin current generated by 
a magnetic moment which precesses in a semi-infinite line of 
lattice sites, calculated at resonance. The surface is at N = 1, 
and the precessing impurity is located at the site TV = 1000. 
In this calculation we set Ef — 0. 



films. If the magnetization in such a film is excited by a 
microwave field, all spins in a given plane parallel to the 
surface and interface will precess in phase. The physical 
situation is thus similar to that in the one dimensional 
model explored in the present paper. We shall explore 
the question of resonant structures in the relaxation rate 
for realisitic three dimensional electronic structures in the 
near future. 

It is our view that the features described in the dis- 
cussion above, while illustrated through use of our one 
dimensional model, are general phenomena that will be 
present in full three dimensional calculations directed 
toward ferromagnetic resonance excitation of coherent, 
spatially uniform spin motions in an ultrathin film on a 
substrate or such a film with a metallic capping layer. 
We note that the existence of quantum well states have 
been directly observed in magnetic multilayers by inverse 
angle-resolved photoemission[28] . 

We now turn to effects specific to our one dimensional 
model. An example of a physical system wherein these 
considerations will be relevant is magnetic moments on 
carbon nanotubes, as discussed in ref. [§1 There is an 
even/odd effect in the spatial distribution of spin current, 
in the presence of a surface. In Fig. |4j with the surface 
site TV = 1, we have illustrated the spatial distribution 
of spin current for an impurity 1000 sites removed from 
the surface, at N = 1001. In Fig. [7J we show what hap- 
pens when the impurity is moved one lattice site closer 
to the surface, /V = 1000. In both figures, the impurity 
is excited at resonance and the calculations assume the 
band is precisely half filled. We see in both figures that 
the spin current vanishes at the surface, as the boundary 
condition in ref. [27] requires. However, when the impu- 
rity is on a site where N is even, there is a discontinuity in 



the spin current at the impurity itself. This is a quantum 
effect; the period of oscillations A for a half filled band is 
two lattice constants, so movement of the impurity one 
lattice constant leads to a phase change of 180 degrees 
for any disturbance it samples from electrons at or very 
close to the Fermi surface. The spin current near the 
impurity on the left side of the figure is a coherent super- 
position of that radiated by the impurity and directed to- 
ward the surface, and that reflected back off the surface. 
There are quantum interference effects associated with 
these two counter propagating currents, so movement of 
the impurity by a single lattice constant affects the total 
spin current samples by it by a large amount. We remind 
the reader that our calculations are all in the limit where 
ballistic transport is envisioned. Such interference effects 
will be absent in a system wherein the electron mean free 
path is shorter than the distance between the impurity 
and the surface. 



IV. CONCLUDING REMARKS 

In this paper we have developed a formalism that al- 
lows the calculation of the spatial dependence of the spin 
current generated by precessing magnetic moments em- 
bedded within or in contact with conducting materials. 
The discussion is carried out entirely within an itinerant 
electron description appropriate to metallic 3d ferromag- 
nets. The magnetic moments whose precession generate 
the spin current are formed through action of on site, 
intra atomic Coulomb interactions of Hubbard character 
as opposed to the localized spin models used in other ap- 
proaches. In this paper, in the interest of simplicity we 
have presented the formalism framed in language suit- 
able to a one dimensional system. The extension to a 
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three dimensional system described by the empirical tight 
binding model we have used in earlier studies [l2 | \T7 \ \18\ 
is very straightforward. We have calculations underway 
directed toward real materials. 

We have also presented a series of numerical calcula- 
tions of the spatial variation of the spin current generated 
by a single precessing magnetic moment embedded in a 
one dimensional lattice of sites, infinite or semi-infinite in 
length. For the case where the moment is embedded in 
an infinite line of sites, we find striking spatial variations 
of the spin current at microscopic length scales that have 
their origin in the Fermi wavelength of the electrons at 
the Fermi surface, and its relationship to the underlying 
lattice constant. These oscillations, as illustrated in Fig. 
[2j are a dynamic analogue of the phenomenon of aliasing 
that enters the description of the static RKKY interac- 
tion between ferromagnetic films separated by a metallic 
film[25]. We can expect such features to be present also 
in three dimensional calculations, it should be remarked. 
For instance, if one considers an ultrathin ferromagnetic 
film on a metallic substrate excited in ferromagnetic res- 
onance, in each layer of the ferromagnet the spins precess 
in phase. In this circumstance, each layer may be viewed 
as a giant spin embedded in a one dimensional lattice. 
The spin current in the substrate will also vary only in 



the direction perpendicular to the interface, and it will 
be constant throughout a given plane. There is thus an 
analogy between this circumstance and the one dimen- 
sional model system explored in the present paper. 

We have also presented studies of the influence of a sur- 
face on the spatial distribution of spin current. Our mi- 
croscopic model produces results in excellent accord with 
the phenomenological discussion of Valet and Fert[S71 
The even/odd effect discussed near the end of Section III 
is fascinating to us. While this is specific to one dimcn 
sional systems, it is our view that in carbon nanotubes 
wherein the spin diffusion length is very long, the consid- 
erations in our discussion are relevant to the description 
of the relaxation rate of a localized moment incorporated 
into such a structure. 



Acknowledgements 

F.S.M.G., A.T.C. and R.B.M. acknowledge support re- 
ceived from CNPq through Instituto do Milenio de Nan- 
otecnologia, Brazil. The research of D.L.M. has been 
supported by the U. S. Department of Energy, through 
grant No. DE-FG03-84ER-45083. 



[1] For a recent experimental and theoretical study, and ref- 
erences to earlier work see J. Lindner, I. Barsukov, C. 
Raeder, C. Hassel, O. Posth, R. Meckenstock, P. Lan- 
deros and D. L. Mills, Phys. Rev. B80, 224421 (2009). 

[2] J. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996); 
ibid. 195, L261 (1999). 

[3] L. Berger, Phys. Rev. B54, 9353 (1996). 

[4] E. Simanek and B. Heinrich, Phys. Rev. B67, 144418 
(2003). 

[5] D. L. Mills, Phys. Rev. B68, 014419 (2003). 
[6] E. Simanek, Phys. Rev. B68, 224403 (2003). 
[7] R. Urban, G. Woltersdorf and B. Heinrich, Phys. 

Rev.Letters 87, 217204 (2001). 
[8] G. Woltersdorf, O. Mosendz, B. Heinrich and C. H. Back, 

Phys. Rev. Letters 99, 246603 (2007). 
[9] F. S. M. Guimaraes, D. F. Kirwan, A. T. Costa, R. B. 

Muniz, D. L. Mills and M. S. Ferreira, Phys. Rev. B81, 

153408 (2010). 

[10] Y. Tserkovnyak, A. Brataas and G. E. W. Bauer, Phys. 

Rev. Letters 88, 117601 (2002). 
[11] M. Zwierzycki, Y. Tserkovnyak, P. J. Kelly, A. Brataas, 

G. E. W. Bauer, Phys. Rev. B71, 064420 (2005); Y. 

Tserkovnyak, A. Brataas, G. E. W. Bauer and B. I. 

Halperin, Rev. Mod. Phys. 77, 1375 (2005). 
[12] A. T. Costa, R. B. Muniz and D. L. Mills, Phys. Rev. 

B73, 054426 (2006). 
[13] K. Lenz, T. Tolinski, J. Lindner, E. Kosubek and K. 

Baberschke, Phys. Rev. B69, 144422 (2004). 
[14] R. Vollmer, M. Etzkorn, P. S. Anil Kumar, H. Ibach and 

J. Kirschner, Phys. Rev. Letters 91, 147201 (2003). 
[15] J. Wang, B.-F. Zhu and R.-B. Liu, Phys. Rev. Letters 

104, 256601 (2010) 



[16] L. K. Werake and H. Zhao, Nature Physics 6, 875878 
(2010). 

[17] A. T. Costa, R. B. Muniz and D. L. Mills, Phys. Rev. 

B70, 054406 (2004). See Fig. 3 for a comparison between 

theory and experiment. 
[18] R. B. Muniz and D. L. Mills, Phys. Rev. B68, 224414 

(2003). 

[19] J. Unguris, R. J. Celotta and D. T. Pierce, Phys. Rev. 
Letters 79, 2734 (1997). 

[20] E. Saitoh, M. Ueda, H. Miyajiima and G. Tartata, Appl. 
Phys. Letters 88, 182509 (2006), H. Y. Inoue, K. Harii, 
K. Ando, K. Sasage and E. Saitoh, J. Appl. Phys. 102, 
083915 (2007), K. Ando, Y. Kajiwara, S. Takahashi, S. 
Maekawa, K. Takemoto, M. Takatsu and E. Saitoh, Phys. 
Rev. B78, 014413, (2008), Y. S. Gui, N. Mecking, X. 
Zhou, G. Williams and C. M. Hu, Phys. Rev. Letters 
98, 107602 (2007) and Phys. Rev. B76, 224430 (2007), 
O. Mosendz, J. E. Pearson, F. Y. Fradin, G. E. W. 
Bauer, S. D. Bader and A. Hoffmann, Phys. Rev. Let- 
ters 104, 046601 (2010), A. Azevedo, L. H. Vilela-Leao, 
R. L. Rodrigues-Suarez, A. F. Lacerda Santos and S. M. 
Rezende, Phys. Rev. B 83, 144402 (2011). 

[21] A. T. Costa, R. B. Muniz, S. Lounis, A. B. Klautau and 
D. L. Mills, Phys. Rev. B82, 014428 (2010). 

[22] A. A. Khajetoorians, S. Lounis, B. Chilian, A. T. Costa, 
L. Zhou, D. L. Mills, J. Wiebe, and R. Wiesendanger, 
Phys. Rev. Lett. 106, 037205 (2011). 

[23] Our notation differs from that used in ref. 1101 and ref. 1111 
These authors use m(t) in place of our n(t). 

[24] R. B. Muniz and D. L. Mills, Phys. Rev. B66, 174417 
(2002). 

[25] D. M. Edwards, J. Mathon, R. B. Muniz and M. S. Phan, 



10 



Phys. Rev. Letters 67, 493 (1991), and J. Phys.: Con- 
densed Matter 3, 4941 (1991); P. Bruno and C. Chappert, 
Phys. Rev. Letters 67, 1602 (1991), and Phys. Rev. B 46, 
261 (1992). 

[26] E. N. Economou, Green's Functions in Quantum Physics, 
Vol. 7 of Springer Series in Solid State Sciences, (Spingcr 
Verlag, Heidelberg, 2006). 

[27] T. Valet and A. Fert, Phys. Rev. B48, 7099 (1993). 



[28] J. E. Ortega and F. J. Himpsel, Phys. Rev. Letters 69, 
844 (1992); J. E. Ortega, F. J. Himpsel, G. J. Mankey and 
R. F. Willis, Phys. Rev. B47, 1540 (1993); F. J. Himpsel, 
J. E. Ortega, G. J. Mankey and R. F. Willis, Advances in 
Physics 47, 511 (1998); J. J. Paggel, T. Miller and T.-C. 
Chiang, Science 283, 1709 (1999). 



